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Abstract 

Identities between Whittaker and modified Bessel functions are derived for partic- 
ular complex orders. Certain polynomials appear in such identities, which satisfy a 
fourth order differential equation (not of hypergeometric type), and they themselves 
can be expressed as particular linear combinations of products of modified Bessel 
and confluent hypergeometric functions. 



1 Introduction 



A class of identities is derived which express Whittaker functions WN,ik(y'^x) 
in terms of modified Bessel functions of the second kind, where k is real, 
is integer or half-integer. In this paper we concentrate on the case where 
A = n + 1/2 where n is a natural number. More explicitly, we will find that 

W^n+i/2,*fc(2x) = xA'^{x)Ky2+ikix) + xA'^*{x)Ky2-ikix) (1) 



where A^(x) is a polynomial of degree n. These polynomials reduce to La- 
guerre polynomials when k = and we will be able to express them as a 
particular linear combination of products of modified Bessel and confluent 
hypergeometric functions. 

We should note that the n = case of this identity was noticed in the solution 
of a "physical" problem; namely the energy eigenfunctions of supersymmetric 
quantum mechanics with an exponential potential [4,5]. 
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2 Proof of identities 



We begin by writing down Whittaker's differential equation [1], 



L{y) = y"(x) + -1 + 



2n + 1 i + 



X 



+ 



y{x) = 0, 



(2) 



which possesses VF„+i/2,tfc(2a;) as a solution. Our strategy will be to show 
that \^{x)Ki/2+ikix) + X'f^* {x)Kif2-ik{x) satisfies this differential equation for 
some polynomial A^(a;) and determine the polynomial as a byproduct. Then 
studying the asymptotics will complete the proof. 

Substituting X'^(x)Ki/2+ik{x) + X'^* {x)Ki/2-ik{x) into the differential equation 
we arrive at 



L{X^Ki/2+ik + c.c.) = 2A.^; K[/2+ik + X^K"i/2+ik + 

I xt — x^ + 



(3) 



-X' + 



X 



x^ 



A„ Ki/2+ik + c.c., 



where +c.c. means add the complex conjugate of the preceding terms. We can 
eliminate the second derivatives i^"i/2+ife(^) by using Bessel's equation, 



K"l/2+ik + -Ky2+ik 



x" 



(4) 



and this gives. 



k' ^'n\j^, , /l/2 + ^/c,fc , l + 2n . fc/A 

An + An Al/2+jfe + C.C. 



X 



X 



2 + 



X 



= L(A^i^i/2+ifc + c.c.). (5) 
Now, we can eliminate the first derivatives -f^i/2+ife(^) using the identities [2] 



xK'^{x) ± vK^{x) = -xK^^x{x)-, 
K,(x) = K_,(x) 



(6) 
(7) 



which lead to 



1 + 2zk^j,f ^ I l + 2ik ^ l + 2n \ 



X 



X^ 



X 



Kl/2+ik + 



(8) 
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The complex conjugate term is not independent. Since -^i/2+ifc(^) = Ki/2-ik{x) 
we can rewrite the whole expression above as {...)Ki/2+ik + c.c. and this can be 
made to vanish if the coefficient of Ki/2+ik{x) is made to vanish; this condition 
corresponds to 



Of course one can consider the complex conjugate version of this differential 
equation, and then we can view them as two linear second-order coupled differ- 
ential equations for and A^*. This will imply that A^ satisfies a fourth-order 
linear ODE. If we substitute X^i^) = J2rrt=o o-m^x"^, we can derive a recurrence 
relation for the coefficients a^^. We find that 



4"^ = 0, (10) 
m(m + 1) (2m - 1) (m + 2ik) (m - 1 - 2iA;)aJ^^2 + (H) 
(1 -I- 2n)m(3m^ + m — 2i/c)a^|i — 

4(1 + 2m){n + m){l + n- m)a^^^ = 0, 1 < m < n - 1. 

This is a rather complicated recurrence relation,; in particular it does not 
generate a hypergeometric series; however given any two members of the se- 
quence it clearly determines the rest. Thus, now we proceed to determine two 
of the coefficients using the asymptotics of the functions. Note that it is actu- 
ally convenient to consider the recurrence relation one gets directly from the 
differential equation above. This is. 



m(m - 2ifc)al;;ii + (1 + '2ri)a^m + (1 - 2m)a^)* = 0, l<m<n (12) 

(n)* (n) p, /-in\ 

and from this one gets to the second order recurrence relation above by elim- 
inating a^^*. It is known [1] that as x — > oo. 



n+l/2 -X 



n 



(14) 
(15) 



which allows us to deduce 
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on+1 

Xl{x) + \l\x) ^ ^x-^\ (16) 



and this tells us that 



providing us with enough information to solve for a^^+i ^-^^d 'we find 

(n) _ ^ 



TT 



Now we turn to the asymptotics for small x. We will need [1,2] 



W^n+i/2,fc(2a;) ^ Jl^!^(2a;)i/2+^'= + c.c. as a; ^ (19) 

Mx) = ^W^,,{2x), (20) 
from which we can derive, 

Xi{x)K,„^.,{x) + C.C ~ 4"^*2'^vf T^^^'^'^'' + ^^^^ 



and upon comparison to (19), obtain 



„<». (22) 



Thus we have derived the first and last coefficient in the polynomial, which 
together with the recurrence relation serve to define A^(x) uniquely. Note that 
the identity is now actually proved as we have shown that A^(x)ii'i/2+jfe(x)+c.c 
satisfies the same differential equation as Wi/2+n,ik{'2x) and possesses the same 
asymptotics and thus they must be the same function. 



3 The polynomials A^{x) 



Let us consider the special case k = for which the identity reduces to a well 
known one. It is clear that in this case the polynomials are real, since both 
and a^^* satisfy the same recurrence relation and boundary conditions. 



4 



Then we see that the polynomials actually satisfy a second order differential 
equation and if we let A°(x) = xyn{2x) we find 

zy\{z) + (1 - z)y'„{z) + ny^iz) = 0, (23) 



and thus we see that yn{z) = CnLn{z) where Ln{z) are the Laguerre polynomi- 
als [2] . Once again, asymptotics can be used to determine the proportionality 
constants c„. Therefore, since Ln{z) ~ 1 as ^ — >• 0, and A°(x) ~ (— l)"'n!x/v^ 
as X ^ 0, we see that c„ = {—l)'^n\/ ^/^^ and thus 

A°(a;) = ^^xL„(2x). (24) 

Now, wc give the fourth order equation that the polynomials satisfy. First 
introduce A^(a;) — xA!^{x); then we have 

xAf + (1 - 2ifc)A^' + (1 + 2n)A^ - 2a;A^'* - Af = 0, (25) 

as our second order equation, and after some work one can eliminate A^* to 
get the rather unsightly answer 

ai{x)A^"" + a2ix)A';'" + a3{x)A'^" + a4{x)A^' + a5{x)A^ = 0, (26) 

a^[x) = x'^[l - 4ik + 4a;(l + 2n)], 

a2{x) = 4x[l - Aik + 3x{l + 2n)], 

a^(x) = -16x^(1 + 2n) + 4x^[l + Aik + Sn{n + 1)] + 

4x(l + 4A;^)(1 + 2n) + 2i{l - 2k){i + k){i + 4A;), 
a4{x) = -32x^(1 + 2n) + Sx[-1 + 2n{n + 1) + 6ik] - 

4(?: + A-)(i + 4A;)(l + 2n), 
a^ix) = 4n(n + l)[4a;(l + 2n) + 3(1 - Aik)]. 

One can work out the indicial equation for this ODE (since there is a regular 
singular point at x = 0) and obtain 

a{a - l)[a^ - a - 4:{1 - k){i + k)] = (27) 

where the solution to the ODE behaves 0. The a = solution of 

course corresponds to our polynomial A^{x). Remarkably, one can write down 
the general solution to this fourth order ODE, which is 

y{x) = CiI^i/2+ikix) Mn+i/2,ik{2x) + C2l-l/2+ik{x) Wn+l/2,ik{2x) 

+C3K_i/2+ifc(a;) Wn+i/2,ik{2x) + CiK_i/2+ik{x) Mn+i/2,ik{2x). (28) 
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Therefore, our polynomial must be such a linear combination and one can use 
the asymptotics as x ^ oo and a; — * to determine all the constants uniquely. 
After a bit of work, one finds the following: 



7rr(l + 2ik) 

= == = i + ^'T(=;rTW' 

2 , , T(—n — ik) ^ 
C3 + -C2 cosh nk + C4 p^_2^^^ = 0' (30) 

TTCs _ r(-^fc)r(i/2 + tk)r{-n + ik) 

' coshTrA; ^V{2ik)V{-n - ik) ' ^ 

which of course can be solved simultaneously. Before doing this let us examine 
the /c = limit. This will lead to C2 = 1, C3 = and C4 = (— l)""''^n!/7r. 
Using [2,3] 



Vr„+i/2,o(2a;) = {-irn\{2xfl^e-^Ll{2x), (32) 

M„+i/2,o(2a:) = {2xf'^e-^Ll{2x), (33) 

/_i/2(x) = W — coshx, (34) 

V T^X 

Ki/2{x) = yJe-^ (35) 



one can show that 



A°(^) = ^-^L:{2x), (36) 
as it should! Finally, solving for C2, C3 and C4 we get. 



ikT{-iky 
2'^i''T{2ik)T{-n-iky 



1 K-'''^) /Q7\ 



2 , , 2iA;r(-iA;)2cosh7rA; 
C3 = -- cosh.^ + ,2..,r(-n-ife)^ + 

r(-iA;)r(-n + iA;) 



V7rr(2iA;)r(l/2 - ik)r{-n - ik) ' 

^ ri-ik)^ri-n + ik) 

^ 27r22^fer(2iA;)r(-n - iA;)^ ^ ' 

which allows us to express the polynomial A^(a;) in terms of modified Bessel 
and confluent hypergeometric functions as follows: 
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^n(^) = C2l-l/2+ik{x) Wn+l/2,jfe(2x) + C3X_i/2+jfc(x) Wn+l/2,ifc(2x) 

+C4i^-i/2+ifc(x) M„+i/2,ife(2a:). (40) 

Now we will give an indication as how to derive the general solution to the 
fourth order equation. Consider the coupled second order equation (25), and 
substitute the trial function A(x) = K_i/2+ik{x)F{x). Using Bcssel's equation 
to eliminate second derivatives of K_i/2+ik{x) and the first order recurrence 
relation (6) to eliminate single derivatives of K_i/2+ik{^), we find that the 
trial function satisfies (25) if 

L{F) = 0, and F*{x) = -F{x). (41) 
Thus possible choices for F{x) are: 

iWn+l/2,ik{^x), 

i(M„+i/2,ifc(2x) + M„+i/2 -ifc(2x)) and 

M„+i/2,ik{2x) - Mn+i/2-ik{2x). 

Hence two solutions to the fourth-order equation are K_i/2+ik{x)Wn+i/2,ik{'^x) 
and K_i/2+ikix)Mn+i/2,ik{'^x)- To get the other two we need to introduce 
I-i/2+ik{x) = I-i/2+ik{x) + Ii/2-ik{x). Notice that this function satisfies the 
same Bessel equation as K_i/2+ik{x) does, and also a very similar first order re- 
currence relation, namely a;/i^/2+ife(^) ~ (~l/2 + ^^)-^-i/2+jfc(a^) = x^-i/2+iki^) 
(this equation has a minus on the RHS for K_i/2+ik{x), see (6)). Substituting 
the trial function I_i/2+ik{x)G{x) leads to a similar condition on G{x) as we 
obtained for F{x): 

L{G) = 0, and G*ix) = G{x). (42) 

Thus /_i/2+jfc(ic)W^„+i/2,ifc(2x) and i^i/2+ik{x)Mn+i/2,ik{'^^) solve the fourth- 
order equation. Using the other two solutions, this means /-i/2+iA;(a;)W^n+i/2,iA;(2a;) 
and I-i/2+ik{x)Mn+i/2,ik{'^x) are also solutions to the fourth order equation, 
and hence we have completed the proof. 



4 Related identities 

We have derived an identity for = n -1-1/2. There also exist similar identities 
for = —n — 1/2. The case = ±n is also of interest; in this case the 
identities look like Wn^iki^x) = y^p„(x)-ft'i+jfc(a;) + ^/xqn{x)Ki}-{x), where 
Pn{x) and qn{x) are polynomials of degree n, but we shall not go through the 
details here. Also note that similar identities probably hold between M^^iki'^x) 
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and /_i/2+ifc(x), since the function Iu{x), like Ki,{x), is symmetric in the order 
and hence /*i/2+jfc(^) = h+{-i/2+ik){x), just hke for K_i/2+ik{x), which was 
an important property we used in the proof. 
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